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We investigate the nonlinear interaction between two photons in a single input pulse at an atomic 
two level nonlinearity. A one dimensional model for the propagation of light to and from the atom is 
used to describe the precise spatiotemporal coherence of the two photon state. It is shown that the 
interaction generates spatiotemporal entanglement in the output state similar to the entanglement 
, observed in parametric downconversion. A method of generating photon pairs from coherent pump 

light using this quantum mechanical four wave mixing process is proposed. 
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^ 1 I. INTRODUCTION 

<: 

! Optical nonlinearities sensitive to individual photons may provide interesting new possibilities of controlling and 
manipulating the quantum states of light 0, 0, IE & IE IE 01 IE 1^- fiol . Possible applications of such nonlinearities 
include quantum nondemolition measurements of photon number pd| and quantum logic circuits for photonic qubits 
12]. Experimentally, sufficiently strong nonlinearities have been achieved in cavity quantum electrodynamics, where 
cavity confinement can enhance the coupling between a single two level atom and the input field 0. By optimizing 
■^j- , the suppression of uncontrollable photon losses in such systems, it may be possible to realize a fully quantum coherent 
t-H ' photon-photon interaction jl.'"i| . The analysis of such a quantum level nonlinearity then requires a quantum mechanical 
t-H i treatment of the spatiotemporal coherence in the input and output fields. Specifically, spontaneous four wave mixing 
effects may entangle the two input photons in their spatial coordinates. This entanglement appears to introduce noise 
in the single photon coherence, even though the two photons are still in a quantum mechanically pure state. 

In order to investigate such effects, we apply a one dimensional model of light field propagation to and from a 
' single two level atom [l4],[l5|. If photon losses are avoided, it is then possible to determine the response functions for 
Oh, single photon and for two photon inputs. Using these response functions, we derive the output state for a resonant 
+-> ' rectangular input. We discuss the implications of this result for coherent input fields and show that it is possible 
to create entangled photon pairs from coherent input light by using an interferometric strategy similar to the one 
^ ■ recently applied in parametric downconversion [l6| . 



II. ONE DIMENSIONAL MODEL OF LIGHT FIELD PROPAGATION 



If the transversal beam profile is known, it is sufficient to describe the propagation of light to and from a system 



using only a single spatial coordinate. In free space, the propagation velocity c is constant. The linear propagation 
process can then be described by a dispersion relation of ui — ck, where k is a scalar |l4j . If this approximation is 
applied to the interaction of electromagnetic field with a single two level system, the transversal profile of the k-space 
eigenmodes is defined by the coupling characteristics of the two level system to the three dimensional field in free 
space. As has been discussed in |14| . the single spatial coordinate r corresponding to the wavevector k then represents 
the distance from the system at r — 0, where negative values indicate propagation towards the system and positive 
values indicate propagation away from the system. 

Figure^ shows a schematic representation of the model. The Hamiltonian of this system can be written as 

Hid = #prop. + #a DS . 

with ^prop. = dk Tick b^bk, 
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FIG. 1: Schematic representation of the one dimensional model for light field propagation in the field-atom interaction. There 
is only one direction of propagation, r < represents light propagating towards the atom and r > represents light propagating 
away from the atom. The interaction takes place locally at r = 0. 



where bk is the photon annihilation operator in fc-space, and cr_ =| G)(E | is the atomic annihilation operator 
describing coherence between the ground state | G) and the excited state | E). The coupling strength is expressed in 
terms of the dipole relaxation rate T. This rate defines the characteristic timescale of the coherent interaction between 
the light field and the two level atom. For convenience, the resonant frequency of the atom has been set to zero. 
Note that this merely corresponds to a rotating frame of reference for the phase oscillations, so that all frequencies 
are expressed as frequency shifts relative to the resonant frequency ujq = cko. 

Experimentally, the model presented here could be realized using a one sided microcavity [l3L ITEf . Losses to 
transversal light field modes can then be minimized and almost all the light emitted by a two level atom inside the 
cavity is emitted along the axis of the cavity. If this ideal condition cannot be met, the model described here could 
still be applied. However, it would be necessary to treat the losses as a transversal mode mismatch between the input 
and output beams and the one dimensional field actually interacting with the single atom. 



It is now possible to formulate the Schrodinger equations for the single photon case by defining the one photon basis 
as | k) for one photon in a fc-eigenstate, and | E) for the excited atom with no photon in free space. The quantum state 
| ip(t)} is then described by the components ip(k;t) = (k \ ip(t)) and i/)(E;t) = (E \ ip(t)). The temporal dynamics of 
these components is given by 



These equations of motion can now be transformed into real space coordinates r by using the Fourier transform 




(1) 



III. LOCAL ABSORPTION AND EMISSION 




(2) 




(3) 



The equation for the propagating field then reads 




(4) 
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As a result of the integration over k, this equation of motion now includes a delta function expressing the locality of 
emission. Since the time evolution should be continuous, this delta function implies a jump of ip( r ]t) at r — 0. By 
integrating equation the discontinuity is found to be given by 

ip(r -> +0;t) - Mr -> -0;f) = -i\ — ip(E;t). (5) 

V c 

Emission and absorption are therefore described by the instantaneous addition of an amplitude proportional to ip(E; t) 
to the single photon wavefunction propagating from r < to r > 0. At r 0, the dynamics of ip(r;t) is simply 
described by linear propagation, ip(r; t) = ip(r — ct; 0). 

In order to obtain the dynamics of ifj(E;t), it is necessary to define the integral corresponding to ip(r — 0;t). The 
proper result is obtained by taking the average of the incoming amplitude ip(r — > — 0;t) and the outgoing amplitude 
ip(r — > +0;i). However, it is convenient to use the result of equation JSJ to express the dynamics of ip(E;t) entirely 
in terms of the incoming amplitude t/j(r — ► — 0;i). It then reads 

j t i>{E- 1) = -T^(E; 1) - iV2cT^(r -> -0; i). (6) 

The amplitude of the excited state ip(E;t) can therefore be obtained from an integration of the previous incoming 
field amplitudes ip(r — ► — 0; t). Since the dynamics of these amplitudes are given by linear propagation at a constant 
velocity of c, they can be obtained from the initial single photon wavefunction at r < using the linear propagation 
dynamics mentioned above. 

With these results, it is possible to integrate the equations of motion from any initial time t- m to any final time 
tout- In particular, the output field within < r < c(t ut — ^in) f° r ^(-^^in) = is given by 

V'Mout) = i } (r-c(t out -t in );t in )-id—'ip(E;t ont -r/c) 
= %l)(r - c{t ont - fin); tj n ) 

-2- [ dr> exp(--(r-r' -c(t out -t in ) ) ^(r';t in ), (7) 



'' • / - c ( t out-*in) V ' 



As the first line of equation Q shows, the output wavefunction is a superposition of a component that propagated 
past the atom unchanged and a component emitted by the excited atom. Since the atom was initially in the ground 
state, the emission can be traced to absorptions of the incoming wavefunction, as represented by the integral in the 
last line of equation JJJ). The output wavefunction at r > can thus be represented as a linear function of the input 
wavefunction at r < 0. 



IV. MANY PHOTON EFFECTS 



The advantage of a local description of the field-atom interaction is that it is easily extended to multiple photons. 
No matter how high the photon density is, we can always define a region from r = — e to r = +e around the atom 
small enough to contain only one photon. In order to solve the field-atom interaction problem for many photons, it 
is therefore only necessary to consider what happens if a photon interacts with the excited atom. 

For this purpose, it is useful to define the many photon Hilbert space as a product space of independent particles. 
The bosonic nature of photons must then be included in the symmetry of the initial state. For reasons of consistency, 
it is then also necessary to distinguish the origin of an excitation, effectively treating the excited state as a state of the 
photon. The two photon wavefunction is then given by the amplitudes for two photons in free space, ^(ri,?^), the 
amplitudes for one photon in free space and one photon at the atom, t/j(ri,E) or ij;(E,r2), and the amplitude for a 
double excitation, tp(E, E). For a two level atom, the latter must always be zero. In the Hamiltonian given by equation 
, this fact is expressed by the difference between the atomic annihilation operator <t_ and the annihilation operators 
of harmonic oscillators. Within the product space of independent particles, this difference is simply represented by 
setting the matrix elements between single excitation and double excitation to zero. The Schrodinger equation for 
the two photon wavefunction then reads 

d Q 

—ip(ri,r 2 ;t) = -c—ip(r 1 ,r 2 ;t) - iV2cT 8(n) ip(E,r 2 ;t) 
at ar\ 

-c—ip{n,r 2 ;t) - iV2cT S(r 2 ) ^(n,E; t) 
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—^(E,r 2 ;t) = -Ti>(E,r2]t)-iV2tfip(ri^-0,r 2 ;t) 

-c^(E,r 2 ;t) (*) 
or 2 

j t i>(n,E;t) = -c-^-ip(n,E;t) (*) 

-r^(r 1 ,£?;t)-iVM ; V(»-i,r2-»-0;t) ) (8) 

where (*) marks the missing two photon absorption terms. This two photon Schrodinger equation describes the nearly 
independent dynamics of two separate photons, except for the absence of absorption for one photon if the other photon 
has been absorbed by the atom. The integration of the two photon Schrodinger equation can therefore be achieved 
by using the single photon results and setting all contributions of double excitation to zero yjj. I n the following, 
however, we will present an alternative solution of the dynamics based on the two photon interaction represented by 
the missing double excitation terms in equation (JSJ . This procedure has the advantage that it can be easily extended 
to three or more photons and may therefore provide a useful foundation for further investigations. 

V. SINGLE PHOTON AND TWO PHOTON RESPONSE FUNCTIONS 

Using the results for local emission and absorption, it is possible to evaluate the effects of the atom-field interaction 
on an arbitrary single photon wavefunction. For this purpose, it is useful to define a time independent characterization 
of the input and output wavefunction. In the context of our model, this characterization is easy to obtain since the 
propagation before and after the interaction processes does not change the shape of the wavepacket. For the single 
photon cases, the input and output wavefunctions can therefore be given by 

V'infa) = , lim = x + ct m^m) 

t in ^-oo 

^out^) = t um , i>{r = x + ct ovA ,;t ou% ). (9) 
t out^+°° 

According to equation Q, the output wavefunction can be obtained from the input wavefunction using a linear 
response function f7i(x;x') such that 

/oo 
dx'U^x'^Jx'). (10) 
-oo 

The single photon response function reads 

f <5(x'-x)-2£exp(-E(x'-x)) for x < x' 

Ui(x;x) = < n t v. / V 11 ) 

l_ for x > x . 

Note that the response function U\(x \ x') is a representation of the unitary operation describing the time evolution of 
the field-atom interaction. It therefore preserves the norm of the wavefunction given by the integral over the absolute 
square. 

Likewise, the field-atom interaction of a two photon wavefunction can be described by a linear response formalism. 
The input and output wavefunctions are then described by 

V> m (xi,x 2 ) = lim ip(n = xi + ct- m ,r 2 — x 2 + ct--t- m ) 
t in -- oo 

^out^i) 2 ^) = t lim, tp(n = X! + ct out ,r 2 = x 2 + ct out ; t out ). (12) 
t out^+°° 

The unitary transform of the input state into the output state can also be described by linear response function, 

/oo 
dx' 1 dx' 2 U 2 {x 1 ,x 2 ;x' 1 ,x 2 )'tpi n (x' 1 ,x 2 ). (13) 
-oo 

If the two photons are always very far apart (xi — x 2 3> T/c), or if the atom is replaced with a harmonic oscillator, 
the propagation of the two photons must be independent of each other. In this case, the response function is equal to 
the product of two single photon response functions, 

{/ lm _(xi,x 2 ;x' 1 ,X2) = C/i(xi;xi)t/i(x 2 ;x 2 ). (14) 
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This response function corresponds to the linear part of the field-atom interaction. However, the absence of two 
photon absorption in the dynamics causes a coupling between the photons that can be described by a nonlinear 
correction AJ7 non ^ m such that 

U 2 (x 1 ,x 2 ;x' 1 ,x' 2 ) = U 1 {x 1 ;x' 1 )U 1 {x 2 ;x' 2 ) + AU n nT1 j - m (xx , x 2 ; x[ ,x' 2 ). (15) 

According to the considerations in the previous section, A[/ non jj n can be found by integrating the contributions from 
double excitations in £^rj n . The result reads 

, J -4^exp(-E(a;i-a:i))exp(-^(x^-a:2)) for Max.{x u x 2 } < Min^, x 2 } 

^nonlin \ Xl > X2 > ^l' x 2> ~ \ l 1 ") 

[ else, 

where the minimum Min{x' l5 x' 2 } effectively defines the latest absorption time and the maximum Max{xi, x 2 } defines 
the earliest emission. Thus the nonlinearity removes all components where the first emission occurs only after the 
second absorption [Til ]. 

It is now possible to derive the output wavefunction for any two photon input wavefunction by integrating equation 
(|T3)l using the expressions for Ui and for A[/ non j- n given by equations (fTTf) and 1 (13) 1 . respectively. If the input is a 
single mode two photon pulse, the input wavefunction can be written as a product state 

tp in (xi,x 2 ) = in (xi)0 in (x 2 ), (17) 

where (j)- m defines the shape of the input pulse. The quantum state of the output field can then be described by 

^out^i^a) =0 out (a;i)'/'out( a; 2) + A^ Ilon ii n .(a;i,a;2), (18) 

where 4> ou i describes the linear single photon response given by 



^out( x ) = / dx'Ux(x;x')(j) in (x'), (19) 
and the nonlinear contribution is directly obtained from 

/oo 
dx[ dx' 2 AC/ nonlirL (xx , x 2 ; x[ , x' 2 ) (f> in (x[ )<t>i n (x' 2 ) . (20) 
-oo 

These equations describe the nonlinear response of the two level atom at the quantum level. It is now possible to 
apply this response function to a variety of input states. In the following, we will focus on the case of a resonant 
rectangular wavepacket. 



VI. THE QUANTUM LEVEL NONLINEARITY AT RESONANCE 

Since the absorption of a photon is strongest at resonance, a resonant input should also produce the strongest 
nonlinear effect in the field-atom interaction. In order to investigate this resonant nonlinearity, we consider the 
response to a rectangular input wavepacket given by 

f -i_ for < x < L 
[ else. 

The linear and nonlinear parts of the output wavefunction for this rectangular wavepacket can be determined analyt- 
ically. They read 



--2=(l-exp(-^))exp(§z) for x < 
r , 



and 



0out(aOH -^(l-2exp(-E(L-a:))) for < x < L (22) 

else 



J -i(l-exp(-^)) exp(^(x 1 + a;2 -2Max{0, a;i ,x 2 })) for x t < L 
A^nonlin.^i'^) = < (23) 
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FIG. 2: Contour plots of (a) the output wavefunction V'out^ 1 ' Xa ) anc ^ 3 ) ^ ne nonlinear component A?/j non j m (a;i,X2) of 
the output for a resonant rectangular input wavepacket of length L — 20c/F. The contour shading corresponds to amplitudes 
ranging from — 4/L for black to +2/L for white. The dark grey shading at the edges of the graphs correspond to zero amplitude. 
The light grey shading of the triangular plateu regions in (a) correspond to an amplitude of 1/L equal to the input amplitude 
of the rectangular wavepacket. 

Figure |21 shows the output wavefunction V'out ana - ^ ne nonlinear component A?/; non i m at an input pulse length of 
L = 20c/r. The most remarkable feature of the nonlinear contribution is its localization around x\ = Xi. This is a 
direct consequence of the local interaction between the two photons. 

A detailed discussion of the two time correlation originating from this spatiotemporal locality of the interaction is 
given elsewhere |l5j . In the present paper, we focus on the coherent properties of the two photon wavefunction. For 
this purpose it is useful to simplify the results by assuming the limit of long pulses, L 3> c/T, and concentrating on 
the region within the pulse, < x% < L. In this limit, the photon-photon interaction becomes independent of the 
pulse shape effects caused by the sudden rise and fall of the rectangular pulse amplitude. The results should then 
apply to any pulse with an input amplitude varying slowly on a scale of c/T, where the pulse length parameter L 
defines the local photon density as 2/L. The output amplitudes are then given by 



1 

L 



(24) 



^nonlin.^i'^) 



4 r 

-exp( \xx 

L c 



X2\) 



(25) 



^outfaii^a) = \ U -4exp(-~|£i - x 2 \) 



(26) 



In the long pulse limit, the linear part of the output wavefunction is nearly equal to the original input pulse. However, 
the nonlinear contribution reduces this overlap by scattering photons into other modes according to 



(V-out I An) = 1+ / dx i dx 2 AW( x i> x 2)A?A non i m (:ri,:E 2 ). 



4 f T 

1 — - / dx~ exp( \x-\) = 1 

L J c 



8C 

TL' 



(27) 



7 



The probability that the two photons will be scattered out of the input mode is therefore approximately equal to 

1 fir 

1-l^outlVWI 2 ^- (28) 

The long pulse limit requires that this fraction is never close to one. However, the result can be used to define a 
scattering cross section for the two photons. If we think of the first photon as being in a random position within the 
pulse, the chance of finding the second photon within a distance < a should be equal to 2a/ L. The interaction cross 
section a for the two photon nonlinearity can then be defined as a = 8c/T. Note that c/T is the coherence length 
of spontaneous emission from the atom. The nonlinear photon-photon interaction mediated by the two level atom 
therefore appears to extend over a region eight times longer than this coherence length. 



VII. ENTANGLEMENT AND FOUR WAVE MIXING IN THE NONLINEAR COMPONENT 

In the long pulse limit, the input mode is very nearly a plane wave resonant with the two level atom (k = 0). It is 
therefore possible to describe the scattering effect as a four wave mixing effect changing the photon frequencies from 
fco = to +k and —k, respectively. The k-space representation of the output wavepacket can be obtained by using 
the local Fourier transform in the spatial region from x% = to Xi = L given by 

^out^i'fe) = y 1 dx 1 dx 2 exp(~ik 1 x 1 )exp(-ik 2 x 2 )ip out (x 1 ,x 2 ) 
L Jo 

~ 4^o4 2 ,0- L (Y2 + C c 2 k 2^ 6 kl,-k 2 , (29) 

where fcj can have values equal to integer multiples of 2tt / L. Note that this discretization of k% is necessary to preserve 
the correct normalization of the quantum state. The phase matching conditions of four wave mixing is expressed in 
equation l|29|) as a Kronecker delta, tffcj,— * 2 , ensuring that the sum of k\ and k 2 is indeed zero. As a result of this strong 
correlation between fcj and k 2l the k-space representation of the two photon output is the Schmidt decomposition of 
the entangled state |17| . 

I Vout) =1 *i = o; k 2 = o) - Yl L(r2 8 + C c2fc2) I fc ; - fc >- (3°) 

According to this representation of the two photon state, the single photon density matrix can be written as a mixture 
of k-eigenstates with 

p = (i - jg) | k = o)(k = o 1 + £ ( £(r2 8 f c2fc2) ) 2 1 *)<* i ■ (3D 

This density matrix defines the single photon coherence of the output. In particular, the frequency spectrum of the 
scattered light is given by a squared Lorentzian, 

Scatter ( fc ) = ~^{ k I P I k ) 



16c 2cT 
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TL 7r(r 2 + c 2 P) 2 



(32) 



Note that the resolution factor Ak = 2n/L is required to adjust the normalization of the continuous spectrum 
^scatter W 1° t ne discrete distribution given by p. Figure |3| shows this scattering spectrum in comparison with the 
spontaneous emission spectrum of the two level atom. It should be noted that the squared Lorentzian of the scattering 
spectrum is narrower than the Lorentzian of spontaneous emission. This spectral feature clearly distinguishes the two 
photon scattering process from an incoherent sequence of absorption and reemission and may serve as an indication 
of spontaneous four wave mixing in experiments where low detection efficiencies prevent an evaluation of two photon 
coincidences. 

As this analysis shows, the resonant nonlinear interaction of the two photons at the atom causes correlated changes 
in the frequencies of the photons. Since the output state is completely quantum coherent, the noise in the single 
photon density matrix actually indicates entanglement between the scattered photons. This situation is quite similar 
to the creation of photon pairs by spontaneous parametric downconversion. In fact, it may also be possible to create 
entangled photon pairs from the spontaneous four wave mixing effect at a single atom nonlinearity by isolating the 
nonlinear part of the two photon response to a coherent input field through destructive interference with an appropriate 
reference pulse. This method will be discussed in the next section. 
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FIG. 3: Frequency spectrum / sca tter °f the photons in the nonlinear component A?/; nomm . The intensity I and the frequency 
k have been scaled in such a way that the area of the spectral line in the graph is equal to one. The dashed line shows the 
Lorentzian line of spontaneous emission from the two level atom derived from the same model Note that the area of this 
line is also one. The comparison shows that the spectrum of photons scatters by spontaneous four wave mixing at the single 
atom is narrower than the spectrum of spontaneous emission. 



VIII. GENERATION OF PHOTON PAIRS USING COHERENT INPUT LIGHT 



Spontaneous four wave mixing can only occur if two photons interact. Moreover, phase matching requires that a 
photon scattered to +k must always be accompanied by a photon scattered to — k. It is therefore possible to use the 
nonlinear photon-photon interaction to generate correlated photon pairs from a coherent input pulse by selecting the 
corresponding output ports in a spectrometer. However, even better results for photon pair creation may be achieved 
if the linear component is removed by interference with another coherent light field using a method similar to the one 
applied to parametric downconversion in |16( . 

For any pulse shape defined by the wavefunction <p, it is possible to define a weak coherent state | a) with a low 
average photon number \a\ 2 <C 1. This coherent state can then be expanded into components with zero, one, and two 
photons. Using | Vac.) for the vacuum state, | 0) for the single photon pulse, and | </>; </>) for the two photon pulse, 
this expansion reads 

|a)«|Vac.)+a|0) + ^|0;0) + ... (33) 

The unitary operator U describing the response of the two level system can now be applied separately to the vacuum, 
to the single photon state, and to the two photon state. The vacuum state is not changed by the interaction at 
all (U | Vac.) =| Vac.)). In the resonant long pulse limit, the single photon component changes its phase by 7r, but 
remains nearly unchanged otherwise. However, the two photon component is changed by the addition of | A^nonlin )■ 
The expansion of the output state therefore reads 

a 2 a 2 
I ^out) ~ I Vac -> - a I 9) + I 4>\ 0) + I A V> n0 nlm.) 

ra 1-") + ^ I A ^nonlin.)- ( 34 ) 

The linear component can therefore be represented by the weak coherent state | —a) with the same coherence 
properties as the original pulse. This coherent pulse can be removed by destructive interference with a much stronger 
reference pulse of the same shape at a high reflectivity beam splitter. Note that the high reflectivity of the beam 
splitter is necessary to avoid quantum noise effects in the interaction that would appear as photon losses in the final 
output. If these conditions are met, the destructive interference may be represented by the displacement operator 
D(a) = exp(a a) — a* a). For |a| <C 1, this operator is only slightly different from 1, but it does have the fundamental 
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property that D(a) \ —a) =\ Vac). The final output therefore reads 

a 2 

D{a) | Vout) » D(a) \ -a) + D{a)-j= | A^ nonlin .) 

« |Vac.) + ^|A^ nonlin > (35) 



This output wavefunction now contains only a zero and a two photon component. The one photon component has 
been eliminated by the interference effects at the high reflectivity beam splitter. It is therefore possible to generate 
entangled photon pairs with a two photon wavefunction described by ^■'4 > Il0Il ii Il using a coherently driven dissipation 
free two level atom and an interferometric setup. The average number of photon pairs created in each pulse is then 
given by 

^T^nonlin. I ^nonlin.) = f^H 4 - ( 36 ) 

In the long pulse limit, it is possible to approximate continuous input light as a sequence of rectangular pulses of 
length L 3> c/T. The intensity of the pump light is then given by I- m = c\a\ 2 / L and the rate of pair creation -Rpgjr is 
given by the average number of pairs per pulse divided by the pulse duration L/c. The result of this estimate reads 

g 

E pair = p 7 hi' ( 3? ) 

where higher order many photon effects are negligible if 7 m <C I\ The pair creation rate is therefore also limited to 
■Rpajj- <C r. However, T is usually in the range of nanoseconds, so considerable pair rates should be possible. 
According to equation (|30|l . the quantum state of the emitted photon pair can be written as 

I A Vwiiin.) = - J2 L(r2 8 + C c 2 fc 2) I fc ; - fc >' ( 38 ) 

where k represents the discretized k-space with Afc = 2it / L. In real space representation, the same entanglement is 
expressed by the coefficients AV> nor j m (xijXa) given by equation These representations show the same time- 

frequency correlations as a phase matched parametric downconversion, that is x\ ~ X2 and k\ = — ki. It may therefore 
be possible to use photon pairs created by four wave mixing in applications similar to those of downconverted photons. 



IX. CONCLUSIONS 



We have described the spatiotemporal dynamics of a one dimensional light field interacting with a single two level 
atom for input states with up to two photons. In the case of a resonant two photon input, the interaction at the atom 
results in spontaneous four wave mixing effects, scattering the photons to higher and lower frequencies. Since this 
scattering effect is fully quantum coherent, the resulting output state is entangled in frequency and time. 

For a coherent state input, it is possible to remove the linear single photon and two photon components by destructive 
interference with a reference pulse. The remaining output then consists of the vacuum state and a small contribution 
from the nonlinear two photon component. This output is very similar to the output from spontaneous parametric 
downconversion. It may therefore be possible to realize a source of entangled photon pairs using the spontaneous four 
wave mixing effects at a single two level atom. 
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